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OPTIMALITY OF NEIGHBOR-BALANCED DESIGNS FOR TOTAL 

EFFECTS 
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The purpose of this paper is to study optimahty of circular neighbor- 
balanced block designs when neighbor effects are present in the model. 
In the literature many optimality results are established for direct ef- 
fects and neighbor effects separately, but few for total effects, that is, 
the sum of direct effect of treatment and relevant neighbor effects. We 
show that circular neighbor-balanced designs are universally optimal 
for total effects among designs with no self neighbor. Then we give 
efficiency factors of these designs, and show some situations where a 
design with self neighbors is preferable to a neighbor-balanced design. 

1. Introduction. In many experiments, especially in agriculture, the re- 
sponse on one plot may be affected by treatments on neighboring plots as 
well as by the treatment applied to that plot. Similarly, in cross-over designs 
the response on one subject in a given period may be affected by the resid- 
ual effects of treatments applied to that subject in the previous periods. In 
both situations optimality results of designs with neighbor balance proper- 
ties are available in the literature. See Shah and Sinha (1989) for a review 
of optimality results of neighbor-balanced cross-over designs with first-order 
residual effects, Hedayat and Afsarinejad (1978), Cheng and Wu (1980), 
Kunert (1984b), Kushner (1997) for cross-over designs, Kunert (1984a) and 
Druilhet (1999) for circular designs. 

However, optimality results are almost always established for treatment 
and neighbor (or residual) effects separately. Usually, one of the aims of 
the experiment is to find a single treatment which can be recommended 
for use on larger spatial areas or over longer time periods than those used 
for individual treatments in the experiment: for example, a single variety of 
wheat to be grown in whole fields, a single drug to alleviate the symptoms 
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of chronic asthma, or a single type of feed to be given to cows throughout a 
whole lactation. When the chosen treatment is in use, its only neighbor will 
be itself; thus the effect of most importance is the sum of the direct effect 
of the treatment and the neighbor effect(s) of the same treatment. 

For cross-over trials, this is the sum of the direct effect and the residual 
effect of the same treatment. In the context of animal feeding trials, Pat- 
terson (1950, 1951) called this the total effect of the treatment; Kempton 
(1991) the permanent effect. Matthews (1988) also used the term total effect 
in cross-over trials, no matter what the area of application. In experiments 
on growing plants, the treatments might be varieties or species. In this case 
the sum of the direct effect of a treatment and its own neighbor effect is 
sometimes called the species effect [McGilchrist (1965)], the variety effect 
[Besag and Kempton (1986)], the monoculture effect [McGilchrist and Tren- 
bath (1971) and Kempton (1982)] or the pure stand effect [Kempton (1985, 
1991, 1997)]. However, the treatments may equally well be pesticides or 
quantities of fertilizer, in which case the French "I'effet plein champ" seems 
more appropriate. However, this terminology is clearly not suitable for drug 
trials. We propose that the sum of the direct effect of a treatment and any 
relevant neighbor effects of that treatment should be called the "total effect" 
irrespective of the area of application. 

In some contexts it has been proposed that the effects of direct and 
neighbor effects are not additive when a treatment has itself as a neighbor 
[Speckel, Vincourt, Azai's and Kobilinsky (1987) for competition between 
sunflowers]. In others [e.g., Kempton (1991)], there are correlations between 
neighboring responses. We do not assume such extra complications here. 

Our designs are in blocks, usually incomplete. Ignoring neighbor effects, 
a block design is inefficient if any treatment occurs more than once in a block, 
so we usually insist that our designs are binary. Thus no treatment is ever a 
neighbor of itself. In particular, it is impossible to achieve orthogonality be- 
tween direct effects and neighbor effects. Some designs optimal separately for 
the estimation of direct treatment and neighbor effects are known in which 
no treatment is preceded by itself. Examples of such designs are given by He- 
dayat and Afsarinejad (1978) for cross-over designs without preperiod, and 
by Magda (1980) for circular cross-over designs. Azais, Bailey and Monod 
(1993) give a catalog of circular neighbor-balanced designs with t — 1 blocks 
of size t or t blocks of size t—1, where t is the number of treatments. These 
latter designs have practical importance in field experiments. Indeed, neigh- 
bor balance implies that the number of replications is divisible by t — 1 and 
experimenters rarely have the resources for 2(t — 1) or more replications. In 
the present paper we aim to show that these designs are universally optimal 
[in Kiefer's (1975) sense] for total effects under models which incorporate 
one-sided or two-sided neighbor effects among the class of all designs with 
no treatment preceded by itself. 



NEIGHBOR-BALANCED DESIGNS 



3 



However, if there are a large number of blocks, then the gain from having 
self neighbors can outweigh the loss from nonbinarity. We use continuous 
block design theory to examine such designs and hence derive efficiency 
factors for binary neighbor-balanced designs. 

2. The designs and the models. All the designs are assumed to be in 
linear blocks, with neighbor effects only in the direction of the blocks (say 
left-neighbor and right-neighbor effects). Because the effect of having no 
treatment differs from the neighbor effect of any treatment, we consider 
only designs with border plots, that is, designs with one plot added at each 
end of each block. The border plots receive treatments but are not used for 
measuring the response variables [see Langton (1990) for the importance of 
such designs] . The plots which are not on the borders are called inner plots. 
The length of a block is the number of its inner plots. We assume that all 
the designs are circular, that is, the treatment on a border plot is the same 
as the treatment on the inner plot at the other end of the block. 

We denote by ^(t,b,k) the set of all circular designs with t treatments and 
b blocks of length k. We assume that k <t and that there is an integer £ 
equal to bk{t{t - 1)}"^ 

Definition 1 . A circular binary block design is a circular design which 
has each treatment at most once in the inner plots of each block. 

Definition 2. A circular neighbor-balanced design (CNBD) is a circular 
binary block design in ^^(^t^h^k) which is a balanced block design in the usual 
sense [Shah and Sinha (1989)] and such that for each ordered pair of distinct 
treatments there exist exactly i inner plots which receive the first chosen 
treatment and which have the second one as right neighbor. 

Definition 3. A circular design neighbor-balanced at distances 1 and 2 
(CNBD2) is a circular neighbor-balanced design such that for each ordered 
pair of distinct treatments, there exist exactly i inner plots that have the 
first chosen treatment as left neighbor and the second one as right neighbor. 

Here are two examples of circular designs neighbor-balanced at distance 
1 and 2. The rows correspond to the blocks and the plots at each end of the 
blocks are the border plots. 
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For a design d, we denote by d{i,j) the treatment assigned to plot j of 
block i: in particular, d{i,0) and d{i,k + 1) are the two treatments applied 
to the border plots of block i. The circularity condition implies that d{i, 0) = 
d{i, k) and d{i, k + l) = d{i, 1). We also denote by Yij the response on plot j 
of block i. All the observations are assumed to be uncorrelated with common 
variance. We deal with two distinct models for the expectation, 

{M 1) E{Yij) = (3i + Td(i,j) + Xd{i,j-i) for 1 < i < & and 1 < j < A:, 

{M2) E(Yij) = (3i + Td{ij) + Ad(ij_i) + Pdii,j+i) 

for 1 <i <b and 1 < j < k. 

The unknown parameters have the following meanings: Pi is the effect of 
block i, T^(ij) is the effect of treatment d{i,j), y^d{i,j-i) is the left-neighbor ef- 
fect of treatment d{i,j — 1), Pd(ij+i) is the right-neighbor effect of treatment 
d{i,j -|- 1). In the standard vector notation, we have 

{Ml) E{Y) = B(3 + TdT + LdX, 

{M2) E{Y) = B/3 + TdT + LdX + RdP, 

where B^Td,Ld and Rd are the incidence matrices of block, treatment, left- 
neighbor and right-neighbor effects. Note that model (A^l) corresponds to 
only one-sided neighbor effect. It is particularly adapted to temporal prob- 
lems with carry-over effects. Model (A^2) corresponds to specific additive 
influence from each left and right neighbor. 

Definition 4. The vector of total effects for models with one-sided 
neighbor effects is defined by (j) = T + \. 

Definition 5. The vector ip of total effects for models with two-sided 
neighbor effects is defined hy if^ = t + X + p. 

3. Some technical tools. We introduce some notation and results used 
throughout the next sections. 

We denote by Ik-, Ik and Jk, respectively, the vector of ones of length /c, 
the (/c, k) identity matrix and the (/c, k) matrix of ones. For any matrix A, we 
denote by the Moore-Penrose inverse of A. The projection matrix onto 
the column span of matrix A is denoted by pr^^). Thus pi'(A) = A{A' A)^ A' . 
We also define pr^^ by pr^^ = / - pr(^). Put Qk = wft^) ^ ~ ^^^-^k- For 
a square matrix A, we denote by tr(j4) the trace of A. For two symmetric 
matrices M and N , M < N means that iV — M is a nonnegative definite 
matrix. A matrix is completely symmetric if it can be written as al + bJ for 
two scalars a and b. 
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Consider the standard partitioned linear model: 

(M) Y = Aa + B(3 + e with E(e) = and Var (e) = a"^!. 

The first lemma is classical [Kunert (1983)], the second, whose proof is in the 
Appendix, gives an upper bound of the information matrix: it generalizes 
a result by Pukelsheim [(1993), page 97] and gives a simple condition for 
equality. 

Lemma 1. Under model (M ) the information matrix C[a.] for the effect 
a is C[a] = A' pr(^) A. 

Lemma 2. Assume that in model (M ) our interest is just for some lin- 
ear combinations of a, say K'a. Then the corresponding information ma- 
trix C[K'a] satisfies C[K'a] < {K'K)+K'C[a\K{K'K)+ with equality if and 
only if C[a] commutes with pr^^). 

Throughout the paper we deal with universal optimality defined by Kiefer 
(1975). A universally optimal design has many good properties: see Shah and 
Sinha (1989) for further details. 

Proposition 3 [Kiefer (1975)]. Assume that a design d* has its infor- 
mation matrix completely symmetric; then d* is universally optimal over a 
class V of designs if and only ifiv{Cd*) = maxdgx)tr(Cd). 

4. Optimality of circular neighbor-balanced designs for total effects. In 

this section, we show that a CNBD (resp. a CNBD2) is universally optimal 
among all the designs with no treatment preceded by itself under the one- 
sided (resp. two-sided) neighbor effect model. 

Lemma 4. Let d* he a circular neighbor-balanced design in ^(^t,b,k) with 
3 < k <t; then the information matrix Cd* \(t>\ for total effects under model 
(Ml) IS Cd* [0] = h{k - 2){2{t - l)r^Qt. 

Lemma 5. Let d* be a circular design neighbor-balanced at distances 
1 and 2 in ft(t,b,k)) with 4 < /c < t. Then the information matrix Cd*[ip] for 
total effects under model {M2) is Cd.[^] = b{k-3) {3{t - l))~^Qt. 

Theorem 6. Under model (Ail ) and for 3 <k <t, a circular neighbor- 
balanced design in i^(^t,b,k) universally optimal for the total effects among 
all the designs with no treatment neighbor of itself. 
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Theorem 7. Under model (A42) and for 4:<k<t, a circular design 
neighbor-balanced at distance 1 and 2 in ^[t,b,k) is universally optimal for 
the total effects among all the designs with no treatment neighbor of itself at 
distance 1 or 2. 

Proof of Lemma 4 and Theorem 6. Put a' = (r'|A'). Then cj) = K'a, 
with K = l2<^ It- 

By Lemma 1 we have Cd[a] = (Td|Lrf)'pr^^(Trf|Lrf) for any design. Here 
K'K = 2It, so, by Lemma 2 we have 
CM < \K'Cd[a]K 

(1) = l{Td prfs) Td + Ta prfj^) + L'^ prf^^ + L'^ VT^fs) ^d} 

= i{4T^ prf^) T, + T^L, + L'^T^ - 2T'M- 
The last equahty comes from the fact that because of the circularity we have 
'^'d^d = L'd^d and B'T^ = E'L^- So pr^^) = pr(5) and so 

Pi'(B) = T'dLd - T'^ pr(5) Ld = T'^Ld - T'^Td + pr^^^ T^. 

Note that (T^Ld)ii is the number of times that treatment i is a neighbor 
of itself at distance 1. Thus, for any design d in ^l(^t,b,k) with no treatment 
neighbor of itself, tr(Crf[(?!)]) depends only on tr(r^pr^-| T^). For a CNBD d* 
in we have T'^,Ld* = L'^,Td* = i{J — I). Hence C^* [a] commutes 

with pr(^) = 2~^ J2 (8> It- Then Lemma 2 gives equality in (1). Because 
T^, pr^Lj) Td* = b{k - l){t - l)~^Qt, (1) establishes Lemma 4 for a CNBD. 
Moreover, a CNBD is a balanced block design, so it maximizes tr(T^ P^'/s) '^'^^ 
among all possible designs of the same size [see Shah and Sinha (1989)] and 
so, by Proposition 3 and Lemma 4, Theorem 6 is established. □ 

Proof of Lemma 5 and Theorem 7. Similarly, under model {M2), 
tp = K'a with a' = {t'\\'\p') and K = It- Because of the circularity, we 
have R'dTd = T'd^d and T'^Rd = L'dTd- Then we have 

(2) Cdm < i{9T>rf5) Td + 2(riLrf + + L'^Rd + R'dLd - &T'M- 
For any design with no treatment neighbor of itself at distance 1 or 2, 

tr(Cd[V']) depends only on tr(T^pr^^ T^). For a CNBD2 d* in ^(^t,b,k)j we 
have R'^tLd* = L'^^Rd* = £{J — I) - Then Lemma 2 gives equality in (2). The 
remainder of the proof is identical to the previous one. □ 

5. Efficiency of circular neighbor-balanced designs. In this section we 
examine efficiency factors of circular neighbor-balanced designs. First we use 
the method developed by Kushner (1997) and Kunert and Martin (2000) to 
construct optimal designs. Even if these designs usually have a large number 
of blocks, they are useful to derive efficiency factors. 
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5.1. Continuous block designs and related upper bounds. For each block u 
of design d, we denote by T^u, L^u and Rdu the incidence matrices cor- 
responding to block u. Thus, = {T'^i \ ■ ■ ■ l^rfb)'! = {L'^^ \ ■ ■ ■ and 
Rd = {R'di \ ■ ■ ■ l-Rrff))'- Inequalities (1) and (2) give, respectively, 

b b 

(3) CM<T.Cdu and Cd[il^]<Y,Cdu. 

u=l u=l 

where Cdu = ji'^T'^uQkTdu + T'd^Ldu + L'd^Tdu - '^T'^yTdu] and 

Cdu = \{^T'^uQkTdu + '^{Tdu^du + L'duTdu) + L'duRdu + R'du^du — ^T^u'^du}- 

Because tic{Cdu) and tr(Cd„) are invariant under permutations of treat- 
ment labels, we may say that two sequences of treatments on a block are 
equivalent if one sequence can be obtained from the other one by relabelling 
the treatments. If we denote by s the equivalence class of the sequence I on 
the block u, we can define 



c(s) = tr(Cd„) 



c(s) = tr{Cdu) 



2 

k 



t t \ 

i=l i=l I 

' 9 * * * \ 

3^ - ;^ zl + 4 E "^^ + 2 E 

V j=l i=\ 1=1 / 



where Ui is the number of occurrences of treatment i in the sequence I, rui is 
the number of times treatment i is on the left-hand side of itself in sequence 
/ and Pi is the number of plots having treatment i on the left-hand side and 
the right-hand side. 

Proposition 8. Consider a design d with h blocks of size k. If s* max- 
imizes c{s), then, under model (A^l), tr(Crf[(/)]) <bc{s*). If s* maximizes 
c{s), then, under model {M.2), tr{Cd[ip]) <bc{s*). 

Proof. Denote by S the number of equivalence classes and by T^dis) 
the proportion of the blocks containing a sequence in the equivalence class 
s. Then we have tr(Crf[(^]) < bJ2s=i'^d{s)c{s) < bc{s*). □ 

5.2. Optimal continuous block designs under model (A^l). Here we char- 
acterize optimal sequences under model (All) and show how to construct 
universally optimal designs without restriction on the competing classes of 
designs. Then we calculate an efficiency factor for a CNBD. 



Notation 1. We denote by [x\ the integer part of the real x and by / 
the real function f{v) = -l + k-v/2-{2-v/k)[k/v\ +v/k[k/v\'^. 
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Proposition 9. Consider model (A^l) and designs with blocks of size k 
(k>3) and t treatments. Then a sequence I* in the class s* that maximizes 
c{s) is characterized by: 

1. The number of different treatments present in the sequence I* is v* , where 
V* maximizes f(v) subject to v £ {2, . . . ,t}. 

2. The number of occurrences of a treatment present in the sequence is either 
n_ = [A;/t;*J or = n_ + 1. 

3. The number of treatments that occur n_ times in I* is v^, where v- = 
k-v* [k/v*\. 

4. The number of treatments that occur n-)_ times in I* is where = 

V* -V-. 

5. Every treatment present in I* has all its occurrences side by side. 
Moreover, c{s*) = f{v*) < k — y/2k with equality if (2A;)^/^ is an integer. 

The proof is given in the Appendix. 

Remark. When k = 2q{q + 1) for some positive integer q, then f{v) is 
maximum at exactly three points: v* = 2q, v* = 2q + 1, v* = 2q + 2. Put 
w = k/[{l + (2k + l)"^/^)/2j. If w is an integer, then f{v) is maximum at 
exactly v = w. If not, then f{v) is maximum at one or both of the two 
integers either side of w. Moreover, it can be shown that, when k is large, 

v*r^[{2ky/^\. 

Example. Table 1 gives the composition of I* depending on k (when 
t>v*). 

For example, for A; = 5 an optimal sequence contains v* = 3 treatments. 
One treatment appears once and two treatments appear twice. So, for in- 
stance, 

/* = (a, b, b, c, c). 

For k = 14: an optimal sequence contains five treatments: one treatment ap- 
pears twice and four treatments appear three times. So an optimal sequence 



NEIGHBOR-BALANCED DESIGNS 

Table 2 



9 



k 3 4 5 
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Eff(d*) 1 1 0.88 
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is, for instance, 

(a, a, b, b, b, c, c, c, d, d, d, e, e, e). 

Note that other optimal sequences can be deduced by circular permutation 
or symmetry. For k = 4 or k = 12, there are three possibilities for /* as seen 
in the remark above. For k = 16 there are two possibilities. 

Theorem 10. Consider designs with b blocks of size k, k > 3, and 
t treatments. Denote by s* an optimal equivalence class of sequences. Then 
a design d* that has each sequence in s* equally often is universally optimal 
among all possible designs with the same size. Moreover, tr(Crf*) = bf{v*) = 
bc{s*) where v* is the number of treatments present in one sequence of s* . 

Proof. By construction, all the sequences in s* are obtained from one 
sequence in s* by relabelling the treatments. Thus, for design d* , Cd* is com- 
pletely symmetric. Moreover, T'^^Ld* = L'^tT^*. So, by the proof of Lemma 
4, the condition for equality in Lemma 2 holds and we have tr (Cd*) = bc{s*). 
So by Proposition 8, d* maximizes the trace and so d* is universally optimal. 
□ 

Consider now the classical criteria <^>p{Cd) = tr(C^ ^/{t-l))^/P, with ti{Mi) = 
tr(M+)"* if g < 0, and $o = hmp_>o ^p- It is weh known that $o, ^i, "^oo cor- 
respond, respectively, to D-, A-, E-optimality [Shah and Sinha (1989)]. More- 
over, $_i(Crf) = l/tr(Cc;). For a completely symmetric matrix C, ^p(C) 
does not depend on p. Thus we can derive efficiency factors for any <I>p of a 
CNBD d* relative to a continuous block design d** constructed in Theorem 
10 by considering Eff(d*) = tr(Cd*)/tr(Cd")> as Table 2 shows. 

When k is large, the efficiency factor for a CNBD in ^(^t,b,k) can be ap- 
proximated by [k — 2){2{k — V2k))~^, which tends to 0.5 when k tends to 
+00. For k = 3 and 4, the efficiency factor of a CNBD is 1, so we have the 
following result: 

Proposition 11. For k = 3 or 4 a CNBD is universally optimal for 
total effects among all possible designs with equal size. 
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Table 3 
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5.3. Efficiency of CNBD2 under model (7W2). In this section we derive 
optimal sequences under model {M2). Unlike the previous section, we can- 
not construct optimal designs from an optimal sequence by considering all 
the treatment relabellings of the initial sequence, essentially because the 
condition for equality in Lemma 2 does not hold. So we just indicate the 
main result without giving the proof and we derive an upper bound for the 
efficiency factors of a CNBD2. 

Notation 2. We denote by / the function 
/>l,^^2) = -1 + A: - 2i;i/3 - 8^2/9 

- (2 - 2vi/k - V2/k) [{k - Vi)/V2\ + V2/k[{k - Vi)/V2\^. 

Proposition 12. A sequence I* in an optimal equivalence class s* is 
characterized by maximizing f{vi,V2), for all possible values attainable where 
vi is the number of treatments appearing once in I* , and V2 is the number 
of treatments appearing at least twice in I* . 

Remark. It can be shown that if I* is an optimal sequence, then not only 
must vi and V2 have values f | and V2 which maximize /, but also all the plots 
receiving the same treatment are placed side-by-side. For A; = 4 an optimal 
sequence contains four different treatments and thus a CNBD2 is universally 
optimal among all possible designs. For /c > 6, an optimal sequence does not 
contain any treatment just once (i.e., f| = 0). As in Section 5.2, we can 
derive the efficiency factor of a CNBD2 (Table 3). 

APPENDIX 

Proof of Lemma 2. Put P = pr^^) = K{K'K)+K' and M = I - P. 
Then a = Pa + Ma and = M. So E{Y) = AK{K' K)+{K'a) + AM'^a + 
B(3. By Proposition 2.3 of Kunert (1983), with his A^, Bi^ and i?2d replaced, 
respectively, by AK{K'K)+, B and AM, we have C[K'a] < {K' K)+ K' A' pvfj^^ AK{K'K)+ = 

{K'K)+K'C[a]K{K'K)+ with equality ifand only if (K'/^)+Er'^'pr(Lj) AM = 0, 
or equivalents {K' K)+ K' C[a]M = 0. If C[a] commutes with P, then C[a] 
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commutes with M, so {K'K)+K'C[a]M = {K'K)+K'MC[a] = because 
K'M = 0. Conversely, if {K'K)+K'C[a]M = 0, then PC[a]M = = MC[a]P, 
so PC[a] = PC[a]{P + M) = (P + M)C[a]P = C[a]P. □ 

Proof of Proposition 9. Let I be a sequence in s and denote by 

V the number of treatments present in If u = 1, then c(s) = 0; thus the 
maximum must be sought on {2,. . . ,t}. Iff > 1, then necessarily rrii <ni — l 
for i such that nj > 0. Thus, X]j ^ ^ ~ ^-^d then 

c{s)<^(^2k-v-^Y.^^=^ (say) 

with equality if nii = rii — 1 for all i such that rij > 1, that is, if all the 
plots containing the same treatment are side-by-side. Fix v. Then, because 
Si=i'^i = k, A is maximum if and only if = [k/v\ ov Ui = [k/v\ + 1 
for any treatment i present in the sequence. So, necessarily, the number of 
treatments present [k/v\ times in the sequence is u[A;/v + Ij — k and the 
number of treatments present [k/v\ +1 times is k — v[k/v\. Points 2-4 of 
the proposition are then established. 

For such a sequence, A = f{v) where the function / is defined in Sec- 
tion 5.2. The function / is continuous (in spite of the integer part). For 
each positive integer p, / is linear on [k/(j>+ l),k/p] with slope (2p^ + 
2p — k)/{2k). The slope increases with p, so / is concave. If there is a posi- 
tive integer q such that k = 2q{q + 1), then the slope is zero on [2^, 2q + 2] 
and any real number in this interval maximizes /. However, v is an inte- 
ger, so f{v) is maximum at exactly three points: v* = 2q, v* = 2q -\- 1 and 
V* = 2q + 2. Otherwise, the slope is never zero, so / is maximized only at 
k/[{l + (2k + l)^/^)/2j = w. If w is an integer, then f{v) is maximum when 

V = w. Otherwise, it is maximum at one or both of the integers on either 
side of w. 

Now, if k/v is an integer, then we have f{v) = {2k — 2k/v — v)/2 = g{v) 
(say). Because g is concave and / is linear on intervals, then for all v, 
f{v) < fl'(^^)- The maximum oi g \s sd, v = (2A;)^/^. Thus, we have f{v*) < 
g{y/2k) = k-V2k. □ 
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